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Isostaticity in two dimensional pile of rigid disks 

Akihiro Kasahara and Hiizu Nakanishi* 

Department of Physics, Kyushu University 33, Fukuoka 812-8581 

We study the static structure of piles made of polydisperse disks in the rigid limit with and 
without friction using molecular dynamic simulations for various elasticities of the disks and pile 
preparation procedures. The coordination numbers are calculated to examine the isostaticity 
of the pile structure. For the frictionless pile, it is demonstrated that the coordination number 
converges to 4 in the rigid limit, which implies that the structure of rigid disk pile is isostatic. 
On the other hand, for the frictional case with the infinite friction constant, the coordination 
number depends on the preparation procedure of the pile, but we find that the structure 
becomes very close to isostatic with the coordination number close to 3 in the rigid limit when 
the pile is formed through the process that tends to make a pile of random configuration. 

KEYWORDS: static structure of granular pile, isostaticity, marginal rigidity, stress propagation 



Isostaticity is a peculiar property that a rigid gran- 
ular pile is expected to have. 1 ) Consider first a pile of 
frictionless deformable grains under gravity. The pile is 
assumed to be stable against small enough external per- 
turbation. Each grain in the pile should be deformed and 
the forces acting through the contacts between the grains 
are elastic force, that is related to the grain deformation 
through elasticity. Imagine that the grains in the pile are 
being hardened gradually. As the grains become harder, 
their deformation diminishes and some of the contacts 
disappear. In the rigid limit, no deformation is allowed 
and the forces acting between the grains should be de- 
termined only by the structure of contact network of the 
pile without any information of deformation. 

Based on such observation, it has been conjectured 1 - 1 
that the contact network in a static pile of non-cohesive 
rigid grains without friction should have a property that 
the network has just enough contacts to be rigid in 
the sense that removal of any contact makes the struc- 
ture flexible. This property is called isostaticity 1 ' 2 ^ or 
marginal rigidity: 3 ' 4 ^ the structure should be unstable 
against an infinitesimal small external force if there are 
less contacts than this, while it should be overconstrained 
with internal stress if there are more contacts. 

It is not obvious, however, whether this holds true 
or not for the pile that is formed through deposition of 
grains that are rigid from the beginning. It is even less 
obvious for the case that the friction between the grains 
comes into the problem. 

A simple way to check whether a pile is isostatic is just 
to count the average coordination number z of grains in 
the pile. Since all the contact forces in the isostatic pile 
have to be determined without redundant constraint, the 
number of independent components of the forces should 
coincide with the number of the conditions for the force 
balance and the torque balance for each grain. This leads 
to z = 2d for a pile of frictionless spheres, z — d(d+l) for 
a pile of frictionless non-spherical grains, and z = d + 1 
for a frictional pile of both spherical and non-spherical 
grains in d dimensions. 



Isostaticity for the piles deposited from a certain ini- 
tial configuration has been examined by numerical simu- 
lations for three dimensional monodisperse spherical sys- 
tems with various elastic constants, 5 ) and it was shown 
that the piles of frictionless sphere becomes isostatic in 
the hard sphere limit, while the pile of frictional spheres 
are more dependent on the way how the piles are pre- 
pared, but are not likely to be isostatic. 

On the other hand, Blumenfeld et al. performed simple 
tabletop experiments for two dimensional non-spherical 
grains of cardboard and concluded that the piles of two 
dimensional grains become isostatic in the limit that the 
number density of a initial configuration and that of a 
final stable configuration coincides. 

In this report, using numerical simulations for two di- 
mensional system, we examine the piles of two dimen- 
sional polydisperse disk system for several ways of prepa- 
ration procedures. 

We performed molecular dynamics simulations on the 
system that consists of two dimensional disks with lin- 
ear elasticity and damping. The system is polydisperse 
with the uniform distribution in the disk diameter with 
the range [0.9, 1] of the maximum diameter ao, and their 
masses are assumed to be proportional to their areas: 





Fig. 1. Two initial configurations: the triangular lattice with 
polydisperse disks (left) and the random configuration with the 
number density 0.6 (right). 
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Fig. 2. The coordination number z for various elastic constant 
k n for frictionless pile (a) and frictional pile (b). The marks rep- 
resent the pile preparation procedure: the Newtonian equation 
with the triangular lattice initial configuration (+), the New- 
tonian equation with the random initial configuration (x), the 
viscous equation with the triangular lattice initial configuration 
(*), the viscous equation with the random initial configuration 
(□). Each mark represents average over six to twelve realizations. 



the mass of the disk with the diameter ao is denoted by 
mo. The bottom of the system is made rough by attach- 
ing the disks with the interval ao, and we employ the 
periodic boundary condition in the horizontal direction. 
The number of the disks N in the system is 400, and 
the horizontal length of the system is 20ao, thus we have 
approximately 20 layers of disks on average. 

Piles are formed by letting the system run under the 
gravitational acceleration g from initial configurations 
until all the disks stop moving. 

The disk at the position x(t) with mass m follows the 
Newtonian equation 



mx(t) = F(t), 



where the force F(t) consists of the gravitation and the 
contact force from the neighboring disks in contact. The 
two particles i and j at Xi (t) and Xj (t) with radii rj and 



Fig. 3. The coordination numbers z for various elastic constant 
k n in the log- log scale frictionless pile (a) and frictional pile (b). 
The same data are plotted using the same marks with those in 
Fig.2 



Tj , respectively, are in contact when the overlap 5ij given 
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is positive. Then the particle j exerts the force Fij 
the particle i 



on 



where and F\j are the normal and tangential com- 
ponents of the force; 



*ij - 



k n 5ijfiij / ~Yn'Vj / j , 



(4) 
(5) 



Here, hij and ty are the normal and tangential unit vec- 
tors, respectively; As^ is the tangential displacement of 
(1) the contact points after the contact, (v\ ■) is rela- 



tive normal (tangential) velocity, k n (k t ) is the normal 
(tangential) elastic constant, and j n ("f t ) is the normal 
(tangential) damping constant. Note that we assume no 
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Newtonian Triangular lattice 
Random config. 


3.97 0.68 

3.98 0.65 


3.15 0.49 
3.09 0.47 


Viscous Triangular lattice 
Random config. 


3.97 0.63 
3.97 0.64 


3.06 0.60 
3.04 0.46 



This consideration prompts us to employ the viscous 
equation 



Table I. The limiting coordination numbers Zoo and the expo- 
nents a for various preparation procedures. 



threshold for the disks to slip during the contact, which 
corresponds to the case with the infinite friction coef- 
ficient. In the case of frictionless disk, we simply set 
F* • = 

In the actual simulations, we use k t — 0.2k n and 7„ = 
[y'mo] for the frictionless case, and k t = 0.2k n and 
In = It — 2\/fcn [\/ m o] for the frictional case. 

As for the initial configurations, we try the triangular 
lattice and the random configuration (Fig.l); the trian- 
gular lattice with the lattice constant <jq is not a regular 
lattice because the disks located at the lattice points are 
polydisperse. The random configurations are prepared by 
putting disks at random with the area fraction around 
0.6. 

Simulations to form piles start from these configu- 
rations with zero particle velocity and stop when the 
kinetic energy of each disks becomes negligibly small, 
namely smaller than 10 _15 [mogcro]- 

Fig. 2 (a) shows the fc„-dependence of the coordination 
number z for the frictionless disks. It can be seen that the 
preparation dependence is very small and z converges to 
the number very close to 4 in the large k n limit for both 
of the initial configurations. The A:„-dependence is well 
represented by the power law 



z - z 



<x k„ 



(6) 



as is shown in Fig. 3(a). 

As for the case of frictional disks, the results are shown 
in Fig. 2(b). There are two things to be noted in compari- 
son with the frictionless case: (i) the discrepancy between 
the two initial conditions shown by + and x is large, (ii) 
the limiting values of the coordination number is sub- 
stantially different from 3, or the value for the isostatic 
structure of the frictional grain in two dimensions. 

The limiting values of the coordination numbers are 
Zoo = 3.15 for the triangular lattice initial condition and 
Zoo = 3.09 for the random initial condition(Tablc I). 

These features indicate that the piles of frictional disks 
are more dependent on the way they are formed than 
those of the frictionless piles. The fact that z^ for the 
random initial conditions is closer to 3 suggests that the 
piles would tend to be isostatic when they are formed 
through the process that retains a random structure. 



7±(t) = F(t) 



(7) 



for the time development during the deposition process 
in stead of the Newtonian equation, in order to see if it 
makes the structure closer to the isostaticity. The viscous 
equation tends to produce random configurations in equi- 
librium because the disks stop as soon as the force bal- 
ance is achieved without the inertia effect. In the simula- 
tion with the viscous equation, we take 7 = 5[m -\/g/(To] 
with 7„ = -ft = in cqs.(4) and (5). 

The data for the piles by the viscous equation are also 
shown in Fig. 3 and Table I. One can see that the limiting 
values for the viscous equation are actually closer to 3 
than those for the Newtonian equation. 

In summary, in this report we have shown the follow- 
ings; The piles of frictionless disks become isostatic when 
the disks are very hard and they are not sensitive to 
the preparation procedure very much, which is consis- 
tent with the conjecture that the pile of rigid grains is 
isostatic. On the other hand, for the piles of frictional 
disks with the infinite friction, the structure depends on 
the preparation process. If the pile is formed from the tri- 
angular lattice with the inertia, the pile structure seems 
to be distinctively different from the isostatic one even in 
the rigid limit, as has been found in the previous work on 
the three dimensional system. 5 ) We found, however, that 
the pile of frictional disks becomes very close to isostatic 
in the rigid grain limit when we employ the deposition 
process that tends to produce random structure. 

The isostaticity of granular pile should play an cru- 
cial role in the stress propagation of granular system. 
In the case of frictionless pile, it has been shown that 
the isostatic structure has some peculiar features: (i) 
unidirectional nature of stress propagation, 6 ) (ii) cor- 
respondence between the force-force response and the 
displacement-displacement response, 2 ' (iii) unstable re- 
sponse that grows as the distance from the point of an 
external perturbation. 1 ) The effects of friction on the me- 
chanical properties of the isostatic structure are now un- 
der investigation. 
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